clc,clear;
x=1800:0.01:2200;
y=9:0.0001:13;
for i=1:40001
Z{1}(i)=Theorem1(1871.2,0.53512,0.30758,1800,x(i));
Z{2}(i)=Theorem1(9.64,0.01643,0.163104,9,y(i));
end
Z{1}=Z{1}';
Z{2}=Z{2}';
% plot(x,Z{1},'r.')
% legend('\mu=0.53512,\sigma=0.30758')
% title('the distribution function of sup_{0<t<Ta}U(t) for Shanghai Composite Index')
% xlabel('sup_{[0,Ta]}U(t)'),ylabel('Probability Density');
plot(y,Z{2},'r.')
legend('\mu=0.01643,\sigma=0.163104')
title('the distribution function of sup_{0<t<Ta}U(t) for Greatoo Intelligent Equipment Inc.')
xlabel('sup_{[0,Ta]}U(t)'),ylabel('Probability Density');